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ABSTRACT 

When  convective  heat  transfer  takes  place  under  large 
temperature  difference  conditions  the  fluid  transport  prop¬ 
erties  may  vary  markedly  over  the  flow  cross-section,  invali¬ 
dating  theoretical  solutions  based  on  constant  properties  and 
empirical  correlations  based  on  small  temperature  difference 
experiments;  This  problem  has  been  investigated  analytically 
for  fully  established  laminar  and  turbulent  flow  in  circular 
tubes,  but  the  only  really  extensive  experimental  work  for  the 
case  of  the  flow  of  a  gas  appears  to  be  for  turbulent  flow 
with  the  gas  being  heated.  The  analytical  solutions  indicate 
that  the  turbulent  flow  heating  problem  shows  the  greatest 
temperature  dependent  properties  effect,  while  the  effects  are 
small  for  laminar  heating  and  negligible  for  laminar  and  tur¬ 
bulent  cooling,  but  experimental  verifies tiorNqf  this  fact  is 
meager  and  some  of  the  simplifying  assumptions  iKthe  analyses 
definitely  reqi^ire^yerif ication.  The  objective  of _the  present 
investigatior^tee^been  to  experimentally  determine  the  effect 
of  large  temperature  differences  for  flow  of  air  in  a  circu¬ 
lar  tube  under  the  conditions  of  laminar  flow  heating,  lami¬ 
nar  flow  cooling,  and  turbulent  flow  cooling. 

Experimental  data  are  presented  for  laminar  flow  heating 
and  cooling  up  to  absolute  temperature  ratios  of  about  1.8, 
and  for  turbulent  flow  cooling  for  absolute  temperature  ratios 
from  1.5  to  2.7.  The  laminar  flow  heating  data,  despite  tem¬ 
perature  rises  of  as  much  as  1100°P,  and  temperature  differ¬ 
ences  as  high  as  550°F,  show  negligible  temperature  dependent 
properties  effects,  and  in  fact  correspond  very  closely  to 
the  Graetz  solution  (extended)  if  all  properties  are  evalua¬ 
ted  at  local  mixed-mean  temperature.  The  laminar  cooling  data 
also  correspond  well  with  the  Graetz  solution  for  mean  Nusselt 
number  with  respect  to  length,  provided  that  properties  are 
evaluated  at  the  logarithmic  mean  temperature  with  respect  to 
end  temperatures  and  wall  temperature.  The  turbulent  cooling 
data  show  no  detectible  effect  of  temperature  ratio  when 
properties  are  evaluated  at  the  logarithmic  mean  temperature. 
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NOMENCLATURE 

English  Letter  Symbols 

p 

A  -  heat  transfer  area,  ft 

A  -  function  of  x+  defined  by  Eq.  (26) 

B  -  function  of  x+  defined  by  Eq.  (26) 

-  constant  in  series  solution,  Eq.  (3) 
c  -  specific  heat  at  constant  pressure,  Btu/(lb°F) 

lr 

D  -  tube  diameter,  ft 

G  -  mass  velocity,  lbs/(hrft2) 

G^  -  function  defined  by  Eq.  (4) 

H' (~T^)  -  function  defined  in  Eq.  (15) 

h  -  convection  heat  transfer  conductance, 

Btu/(hrft2oF) 

I  -  electrical  current,  amps 

k  -  thermal  conductivity,  Btu/(hrf t2oF/f t ) 

L  -  tube  length,  ft 

L  -  dummy  nondimensional  tube  length  variable  of 
same  character  as  x+ 

o 

L  -  heat  leak,  Btu/in 

P  -  power  input,  Btu/in 

q  -  heat  flux,  Btu/(hrft2) 

q1  -  heat  transfer  rate  per  unit  length,  Btu/(hrft) 

p 

q"  -  heat  flux  (same  as  q  ),  3tu/(hrft  ) 

R  -  electrical  resistance,  ohms 

Rn  -  eigenfunctions,  see  Eq.  (3) 

r  -  radial  coordinate  measured  from  tube  centerline ,ft 

r  -  tube  radius,  ft 

o 

r+  -  nondimensional  radial  position,  r/r 

o 
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T 

T 


temperature,  R 


o 


uniform  fluid  temperature  at  beginning  of  tube,  R 

oT 


w 


T 


m 


t 

U 

V 

w 

X 


wall  surface  temperature,  R 
mixed-mean  fluid  temperature,  °R 
thickness  of  tube  wall,  ft 

overall  heat  transfer  conductance,  Btu/(hrft°F) 
fluid  velocity,  ft/sec 
fluid  flow  rate,  lbs/hr 
axial  coordinate  in  tube,  ft 


nondimensional  axial  coordinate,  (x/rQ)/(N 


ReNPr 
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Greek  Letter  Symbols 


a 

a 


thermal  diffusivity,  ft/hr£ 


temperature  coefficient  of  electrical  resistivity, 
l/°F 


r 


m 


€ 

€H 

0 

0„ 


m 


'n 


-  eigenvalue  defined  by  Eq.  (l4) 
heat  exchanger  effectiveness 
turbulent  eddy  diffusivity,  ft/hr2 
nondimensiona.l  temperature  defined  by  Eq.  (2) 
nondimensional  mixed-mean  temperature 
eigenvalue  defined  by  Eq.  (5) 

-  coefficient  of  viscosity,  lbs/(hrft) 


Nondimensional  Groupings 


N 


Nu 


N. 


Nu, 


-  Nusselt  number,  hD/k 
local  Nusselt  number 


N. 


Nu 


m 


mean  Nusselt  number  with  respect  to  tube  length 


%u  “  Nusselt  number  based  on  constant  fluid  properties, 
iso  i.e.,  with  temperature  difference  approaching 


zero . 

N^  -  Reynolds  number,  DG/p. 
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-  Prandtl  number,  /k 

-  number  of  beat  transfer  units  in  a  heat 
exchanger,  AU/(wCp) 

-  nondimensional  tube  length  parameter, 

(x/ro)/(NReNpr) 
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SECTION  I 

INTRODUCTION  AND  OBJECTIVES 


The  question  of  the  influence  of  temperature-dependent 
property  variations  on  the  heat  transfer  characteristics  of 
a  gas  flowing  in  a  circular  tube  arises  whenever  the  tempera¬ 
ture  difference  for  convection  becomes  more  than  100-200°F. 
Since  most  comection  heat  transfer  data  have  been  taken  from 
experiments  involving  small  temperature  differences  (a  temp¬ 
erature  difference  being  small  if  the  mass  and  thermal  trans¬ 
port  properties,  such  as  viscosity  and  thermal  conductivity, 
do  not  vary  significantly  through  the  flow),  there  is  little 
ambiguity  involved  when  this  data  is  correlated  by  the  use 
of  dimensionless  parameters  Involving  the  mass  and  thermal 
transport  properties  of  the  fluid.  However,  the  following 
question  arises  when  applying  these  data  to  situations  in¬ 
volving  large  temperature  differences:  at  what  temperature 
should  the  properties  be  evaluated;  and  what  influence  has 
a  large  property  variation  across  the  tube? 

In  the  case  of  small  temperature  differences  this  ques¬ 
tion  does  not  arise  because  the  transport  properties  are  es¬ 
sentially  constant,  both  across  the  flow  cross-section  and 
in  the  direction  of  flow. 

For  large  temperature  differences,  it  is  usually  most 
convenient  to  evaluate  the  local  properties  at  the  local 
mixed-mean  temperature.  But  if  the  surface  conductance  is 
evaluated  from  experimental  correlations  constructed  with 
small  temperature  difference  data,  and  if  all  transport  prop¬ 
erties  are  evaluated  at  the  local  mixed-mean  temperature, 
will  the  actual  local  surface  conductance,  in  the  case  of 
large  temperature  difference,  be  higher  or  lower,  and  if  so 
how  much? 

It  should  be  noted  here  that  temperature  dependent  prop¬ 
erty  variation  across  the  flow  cross-section  is  not  the  only 
Important  variation.  Essentially  the  same  questions  as  were 
raised  about  variations  across  the  flow  cross-section  can  be 
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raised  ab  -t  variations  In  the  direction  of  the  flow.  In 
fact,  frequently  the  two  effects  cannot  be  separated,  as  In 
the  case  of  Intermediate  effectiveness  heat  exchangers.  It 
Is  to  be  expected  that  for  low  effectiveness  heat  exchangers 
the  variation  across  the  flow  cross-section  will  be  most  Im¬ 
portant,  and  that  for  high  effectiveness  heat  exchangers  the 
variation  lr  the  flow  direction  will  be  more  significant. 

As  the  temperature  differences  encountered  In  applica¬ 
tions  have  increased,  this  question  of  property  variation  has 
become  more  significant.  There  is  a  moderate  amount  of ^ex¬ 
perimental  data  involving  large  temperature  differences  avail¬ 
able,  and  several  methods  of  correlation  have  been  tried  with 
varying  degrees  of  success.  A  common  approach  is  to  evaluate 
certain  or  all  of  the  properties  which  occur  in  the  dimension¬ 
less  parameters  at  some  temperature  intermediate  between  the 
mixed-mean  temperature  and  the  tube  wall  temperature,  and  some 
at  the  mixed-mean  temperature.  In  the  case  of  gases  this  may 
involve  abandonment  of  the  mass  flow  velocity  (G),  a  useful 
parameter  about  which  there  is  no  ambiguity. 

In  the  case  of  convection  heat  transfer  to  or  from  liq¬ 
uids,  density  and  thermal  conductivity  can  usually  be  treated 
as  constants,  as  the  variation  of  these  properties  is  usually 
small.  It  is  thus  only  the  viscosity  which  varies  markedly, 
and  hence  the  effects  of  temperature  dependent  property  varia¬ 
tion  on  the  heat  transfer  characteristics  are  simply  the  ef¬ 
fects  of  a  deformation  of  the  velocity  profile.  In  this  case 
a  correlation  in  which  the  viscosity  is  evaluated  at  an  in¬ 
termediate  or  "film"  temperature,  or  the  use  of  a  wall- 
vlscosity-to-mixed-mean-viscosity  ratio  as  a  correction  factor 
is  sufficient. 

In  the  case  of  convection  heat  transfer  to  gases  the  sit¬ 
uation  is  more  complex.  Not  only  can  viscosity  vary  markedly 
across  the  flow  cross-section  and  in  the  flow  direction,  but 
also  the  density  and  the  thermal  conductivity.  In  effect, 
with  gases,  the  momentum  and  the  energy  equations  are  coupled 
and  must  be  solved  simultaneously.  It  is  no  longer  possible 
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to  easily  predict  the  effects  of  the  variation  of  the  tempera¬ 
ture  dependent  transport  properties.  As  an  illustration  of 

this  point,  and  as  an  illustration  of  the  sensitivity  to  ge- 

1* 

ometry,  the  analytical  solutions  which  do  exist'  indicate 
that  for  laminar  flow  between  parallel  plates  the  surface 
conductance  will  increase  when  large  temperature  differences 
are  involved  in  the  case  of  heating,  while  for  laminar  flow 
in  circular  tubes  under  similar  conditions  the  surface  con¬ 
ductance  will  decrease.  Fortunately,  the  important  trans¬ 
port.  properties  of  most  gases  vary  in  a  similar  fashion  with 
temperature,  and  hence  it  is  to  be  expected  that  experimental 
data  obtained  with  air  as  the  fluid  medium  would  apply  to 
other  gases  as  well. 

Most  of  the  experimental  data  which  involves  large  tern- 

O 

perature  differences  with  gases  has  been  taken  from  experi¬ 
ments  in  which  air  in  turbulent  flow  was  heated .  This  data 
can  be  successfully  correlated  by  a  correlation  of  the  form, 

NNu  ■  ^NRe-NPr>(^)a 

where  a  is  an  empirical  constant,  and  where  all  properties 
are  evaluated  at  the  mixed-mean  temperature  of  the  fluid. 

The  areas  of  greatest  need  now  for  experimental  data  are  tur¬ 
bulent  flow  with  the  gas  being  cooled,  and  laminar  flow,  both 
heating  and  cooling. 

The  experimental  data  taken  in  this  work  has  been  carried 
out  to  help  fill  this  gap.  In  particular,  the  objectives  of 
this  work  include  the  following: 

1.  To  experimentally  determine  the  effect  of  large  tem¬ 
perature  differences  for  a  gas  being  heated  at  approximately 
constant  heat  rate  per  unit  of  tube  length  in  laminar  flow 
in  a  circular  tube,  and  to  determine  if  an  entry-length  solu¬ 
tion  of  the  Graetz  type  (based  on  constant  properties)  can  be 
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used,  with  suitable  modifications  if  necessary,  under  largp 
temperature  difference  conditions.  This  part  of  the  inves¬ 
tigation  is  concerned  with  local  heat  transfer  conductances 
measured  along  a  tube. 

2.  To  experimentally  determine  the  effect  of  large 
temperature  differences  for  a  gas  being  cooled  with  a  con¬ 
stant  surface  temperature  in  laminar  flow  in  a  circular  tube. 
This  part  of  the  investigation  is  concerned  with  mean  heat 
transfer  conductances  with  respect  to  tube  length. 

3.  To  experimentally  determine  the  effect  of  large  tem¬ 
perature  differences  for  a  gas  being  cooled  with  a  constant 
surface  temperature  in  turbulent  flow  in  a  circular  tube. 

This  part  of  the  Investigation  is  concerned  with  mean  heat 
transfer  conductances  with  respect  to  tube  length. 

The  following  will  be  divided  into  six  sections.  Be¬ 
cause  of  the  great  utility  of  the  available  analytical  solu¬ 
tions,  and  because  they  are  used  as  a  theoretical  reference 
with  which  to  compare  the  experimental  data,  the  second  sec¬ 
tion  summarizes  these  solutions.  The  third  summarizes  and 
discusses  previous  investigations  of  the  problem  of  tempera¬ 
ture  dependent  property  variations  in  gases  flowing  in  circu¬ 
lar  tubes.  The  fourth  section  is  a  description  of  the  experi¬ 
mental  apparatus.  The  fifth  summarizes  the  data  reduction 
procedure.  In  the  sixth  the  experimental  results  are  dis¬ 
cussed,  and  the  conclusions  are  summarized  in  the  seventh. 


SECTION  II 

CONSTANT  PROPERTY  ANALYTICAL  SOLUTIONS 


In  this  section  a  brief  summary  of  the  available  analyti 
cal  solutions  for  heat  transfer  to  a  fluid  flowing  in  a  cir¬ 
cular  tube  is  presented.  More  detailed  information  can  be 
found  in  the  given  references.  All  the  solutions  presented 
are  for  the  case  of  constant  properties  and  fully  developed 
velocity  profile.  Thus  they  are  only  expected  to  apply  -where 
the  temperature  differences  involved  are  small,  and  where  the 
velocity  profile  is  fully  developed  at  the  start  of  the  heat¬ 
ing  section,  or  where  the  Prandtl  number  is  large  so  that  the 
velocity  profile  develops  very  much  more  rapidly  than  the 
temperature  profile. 

The  Graetz  solution  covers  the  case  of  fully  established 
laminar  flow  (parabolic  velocity  profile  throughout)  with  a 
uniform  fluid  temperature  at  the  point  where  the  wall  surface 
temperature  suddenly  steps  from  the  fluid  temperature  to  a 
new  temperature  and  remains  constant  thereafter.  The  solu¬ 
tion  then  yields  the  fluid  temperature  as  a  function  of  a 
radial  position  r+  ,  and  an  axial  position,  x-1"  .  These 
results  are  obtained  starting  with  the  differential  energy 
equation . 


Id  fr,  dT-,  _  V  dT 
r  3F  lr  3FJ  ~  a 


(1) 


After  Introduction  of  the  parabolic  velocity  distribution, 
Eq.  (l)  may  be  rendered  nondimensional  as  follows: 


1 _ d_ 
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+  de 


dr 


/ ,  +2  \ 
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The  boundary  conditions  for  (2)  are 

6(0,r+)  1 

e(x+,i)  =  0 

Equation  (2)  with  the  given  boundary  conditions  is  an  eigen¬ 
value  problem,  and  the  solution  may  be  written  in  the  form 


00 
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With  the  temperature  distribution  established  the  heat 

from  T31  P03ltl°n’  q(X+)  '  may  be  -adily  computed 

irom  ohe  gradient  at  tne  wall. 


where 
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Llpkis1  has  computed  a  large  number  of  the  eigenvalues 
and  constants  for  Eq.  (<,).  Sellar5,  Trlbus<  ^ 

Obtained  an  asymptotic  solution  which  is  accurate  for  large 

n.  The  asymtotic  expressions  for  the  eigenvalues  and  con 
starts  are, 

*n  =  Hn  +  8/3  (5) 

Cn  "  1-012T6A-V3  (6) 

at  flUl<3  temperature’  em  -  can  be  evaluated 

any  x  by  integrating  Eq.  (h)  with  respect  to  x+  and 

ng  an  energy  balance  on  the  flow  stream.  The  result  is 
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at  allTiT  teat  flUX  3nd  mlXSd-raean  CemPerature  established 
Points,  a  local  heat  transfer  coefficient  and  a  local 
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Nusselt  number  can  readily  be  defined  and  evaluated. 
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Zor  large  x  the  series  are  rapidly  convergent,  and  even 
tually  only  the  terms  corresponding  to  n  =  0  are  signifi¬ 
cant.  Then  Eq.  (8)  reduces  to  the  expression  for  the  asymp¬ 
totic  Nusselt  number  for  a  long  tube. 
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Using  the  Graetz  step  function  solution  as  a  base,  the 
method  of  superposition  can  now  be  employed''1  to  solve  two 
types  of  variable  surface  temperature  problems:  (a)  the  sur¬ 
face  temperature  distribution  with  length  is  prescribed  and 
it  is  desired  to  evaluate  heat  flux  at  all  points;  (b)  the 
heat  flux  distribution  with  length  is  prescribed  and  it  is 
desired  to  evaluate  the  surface  temperature  at  all  points. 

If  the  surface  temperature  variation  Is  prescribed,  the 
heat  flux, 

4  k 

q(*  )  =7 


q(x+)  ,  at  any  point  may  be  evaluated  from 
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er+(x+-L,i)ATw 

(10) 

(11) 


where 


-*‘(x+-L) 


V  =  -2  S 
n=0 


Gne 
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L  is  a  dummy  variable  varying  from  0  to  x’  .  The  inte¬ 
gral  in  Eq.  (10)  represents  the  effects  of  a  continuous  vari¬ 
ation  in  the  surface  temperature  while  the  summation  accounts 
for  steps  or  discontinuities.  Thus  ATWl  is  the  height  of 
the  1th  step  located  at  . 

If  the  heat  flux  is  prescribed  and  the  wall  temperature 
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Is  to  be  solved  for,  the  problem  is  more  complicated.  Only 
the  computing  equations  v/ill  be  given  here.  For  more  de¬ 
tailed  information,  see  reference  3  .  The  computing  equa¬ 

tions  are  as  follows: 

T  -  T  =  £2.  /  g(x+  -  L,l)q(L)dL  (12) 

w  O  K  Q 


wnere  q(L)  is  the  prescribed  heat  flux  as  a  function  of 
the  dummy  variable  L  ,  and  the  function  g(x+,l)  is  given 
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The  values  of 
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y  are  the  zeros  of  the  equation. 
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where  s  is  the  Laplace  transform  variable.  H'(-ym)  is 

then  computed  after  differentiating  Eq .  (l4)  and  substitut- 
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ing  values  of  y 
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In  reference  3  the  first  three  terms  for  evaluation  of  Eq . 

(lp)  have  been  computed  and  tabulated.  More  recently  Sparrow 
h 

and  Hallman  have  solved  this  problem  directly  rather  than 
through  a  transform  of  the  constant  wall  temperature  problem. 
They  have  evaluated  the  first  five  eigenvalues  (y  )  ana  have 
obtained  an  asymptotic  solution  for  yr  . 

c  >fl 

Sleicher  and  Tribus  treat  the  turbulent  flow  counter¬ 
part  of  the  Graetz  laminar  flow  solution,  i.e.,  fully  es¬ 
tablished  turbulent  flow  with  a  step  change  in  surface  tem¬ 
perature.  The  differential  equation  solved  is  the  same  energy 
equation,  Eq.  (l),  used  for  the  laminar  flow  problem,  the 
boundary  conditions  are  Identical,  and  the  solution  car  be 
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put  into  the  same  form  and  nomenclature.  The  problem  is  dif¬ 
ferent  from  the  laminar  flow  problem  in  that  the  velocity  ’ n 
Eq.  (l)  is  now  the  turbulent  velocity  profile  instead  of  para¬ 
bolic,  and  the  thermal  diffusivity,  a  ,  is  replaced  by 
(a  +  eu)  ,  where  €„  is  the  turbulent  eddy  diffusivity  for 
heat  transfer.  The  turbulent  velocity  profile  varies  with 
Reynolds  number,  and  the  eddy  diffusivity  is  a  function  of 
radial  position,  Reynolds  number,  and  Prandtl  number.  The 
problem  is  thus  more  complex,  and  the  eigenvalues  and  con¬ 
stants  will  be  functions  of  the  Reynolds  number  and  Prandtl 
number.  Sleicher  and  Tribus  used  empirical  data  for  the  ve¬ 
locity  profile  and  an  empirically  based  modification  of  the 
Jenkins  analysis  for  the  eddy  diffusivity.  Their  results 
can  be  used  to  solve  the  turbulent  flow  problem  where  the 
surface  temperature  is  any  prescribed  function  of  tube  length 
through  the  use  of  Eqs.  (10)  and  (11 ). 

The  Sleicher  and  Tritus  results  3rc  limited  in  two  re¬ 
spects.  First,  they  are  limited  to  Pranuil  numbers  equal  to 
or  less  than  unity.  From  the  point  of  view  of  the  prescribed 
surface  temperature  problem,  this  is  not  serious,  because  only 
at  low  Prandtl  numbers  does  a  variation  in  surface  temperature 
markedly  affect  the  surface  conductance.  The  second  limita¬ 
tion  is  that  only  the  first  three  eigenvalues  and  constants 
have  been  evaluated.  In  some  technical  problems  these  are 
insufficient  to  evaluate  the  series  which  arise.  Also,  it  is 
not  possible  to  handle  the  technically  important  problem  where 
the  heat  flux  is  prescribed  without  higher  order  eigenvalues 
and  constants.  Sparrow,  Hallman,  and  Siegelc  have  solved 
Eq.  (l)  for  turbulent  flow  in  circular  tubes  for  Prandtl  num¬ 
bers  of  0.7,  10,  and  100  for  the  case  of  constant  heat  rate 
per  unit  of  tube  length. 

Q 

Kays  and  Nicoll  demonstrated  that  a  simple  extrapola¬ 
tion  of  tne  Sleicher  and  Tribus  results  could  be  used  to  es- 

/  2 

tablish  the  higher  order  terms  (KJ:  and  G  )  with  an  accuracy 
sufficient  to  evaluate  the  yu  and  H'(-y2)  to  an  accuracy 
commensurate  with  the  uncertainty  involved  in  the  knowledge  of 
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the  turbulent  eddy  diffusivity  at  low  Prandtl  numbers. 

The  results  of  all  these  analyses  have  been  summarized 

in  Appendix  A.  These  solutions  are  used  as  the  theoretical 

reference  to  which  to  compare  the  experimental  data  collected 
in  tnls  work. 


SECTION  III 

SUMMARY  OF  PREVIOUS  WORK  ON  THE  INFLUENCE  OF 
TEMPERATURE  DEPENDENT  GAS  PROPERTIES 


A  search  of  the  literature  indicates  that  only  two  at¬ 
tempts  have  teen  made  10  analytically  predict  the  effects  of 
temperature  dependent  property  variations  on  the  heat  trans¬ 
fer  characteristics  for  gas  flowing  in  circular  tubes.  These 

10  7 

are  tne  work  of  Diessler  and  the  work  of  Sze  .  Both  give 
essentially  the  same  results.  Botn  authors  used  essentially 
the  same  simplifying  assumptions.  The  basic  differences  be¬ 
tween  the  two  are  ir.  the  assumptions  made  atout  tb/.  functional 
relationship  between  the  properties  and  temperature.  Diessler 
assumes  that  both  the  thermal  conductivity  and  the  viscosity 
vary  witn  tne  absolute  temperature  to  the  0.6B  power. 
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and  that  the  specific  heat  is  constant,  i.e.,  the  Prandtl 
number  is  constant.  Sze  used  the  experimental  variation  of 
these  properties  with  temperature  for  air,  and  thus  it  is 
felt  that  his  results  are  more  pertinent  to  the  experimental 

vi  a  i  i  1  4-'-'  r  ^  •  -v  1  > 

O  a.  wux  J  n  U  i  /\  # 

Both  Sze  and  Diessler  assumed  fully  established  velocity 
and  temperature  profiles  and  constant  heat  rate  per  unit 
length.  It  is  to  be  expected  that  the  effects  of  property 
variation,  like  the  surface  conductance,  will  depend  at  least 
to  some  extent  upon  the  variation  of  heat  flux  with  length, 
although  this  may  well  be  a  second  order  effect. 

The  calculating:  procedure  used  for  laminar  flow  was  an 
Iterative  procedure  in  which  a  velocity  profile  was  assumed 
and  then  a  temperature  profile  calculated.  Then,  using  this 
temperature  profile,  the  variation  of  viscosity  and  thermal 
conductivity  with  radial  position  was  determined.  With  the 
variation  of  the  viscosity  with  radius  a  new  velocity  pro¬ 
fil‘d  was  calculated,  and  with  this  velocity  profile,  ana  with 
the  variation  of  thermal  conductivity  with  radial  position,  a 
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new  temperature  profile  was  calculated.  This  process  was  re¬ 
peated  until  the  results  were  unchanged. 

From  the  temperature  profile  and  the  velocity  profile  the 
mixed-mean  temperature  can  be  determined.  Then,  using  the  de¬ 
fining  equation  for  the  surface  conductance,  anu  evaluating 
the  thermal  conductivity  at  the  tnixed-mean  temperature,  a 
Nusselt  number  can  be  determined. 

For  turbulent  flow  essentially  the  same  method  was  used, 
although  assumptions  had  to  be  made  about  the  variation  of  the 
eddy  dif f usivitles  as  a  function  of  radial  position. 

Sze's  results,  plotted  on  Figs,  la  and  lb,  can  be  approxi¬ 
mately  represented  or.  logarithmic  plots  by  straight  lines,  in¬ 
dicating  that  an  experimental  correlation  of  the  form. 
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Nu  Nu 
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should  suffice.  The  idea  of  a  straight-line  correlation  is 
actually  better  substantiated  by  the  only  really  significant 
experimental  data,  also  shown  on  Fig.  12.  It  will  be  noted 
that  for  turbulent  heating  the  NACA  experimental  data  lie  very 
close  to  a  straight  line  with  a  =  -0.575.  The  analysis  lies 
very  close  to  the  experimental  data  for  T  /T  about  1.5.  but 
indicates  a  less  pronounced  effect  for  temperature  ratios  be¬ 
tween  2  and  3*  The  analysis  involved  an  iterative  procedure 
which  lost  in  precision  at  higher  temperature  ratios,  and  this 
is  believed  to  possibly  be  the  reason  for  the  higher  tempera¬ 
ture  ratio  discrepancy.  The  experimental  accuracy  of  the  NACA 
data  is  probably  not  sufficiently  high  to  warrant  three  places 
in  the  exponent,  but  a  =  -0.5  corresponds  reasonably  closely 
to  both  analysis  and  experiment  and  is  recommended  as  a  good 
approximation  for  turbulent  heating. 

For  turbulent  cooling  the  effect  is  evidently  consider¬ 
ably  less  than  for  heating,  and  this  is  also  true  for  laminar 
flow,  and  the  available  laminar  boundary  layer  solutions  in¬ 
dicate  the  same  trend.  A  line  for  a  =  -0.15  has  been  drawn 
on  Fig.  la,  but  since  there  are  only  two  points  (other  than 
the  origin)  upon  which  to  base  a  line,  a  smaller  exponent 
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could  be  Just  as  well  Justified. 

For  laminar  flow  there  are  no  experimental  data  for  com¬ 
parison.  Sze 1 s  results.  Fig.  lb,  would  indicate  a  moderate 
decrease  in  conductance  (or  Nusselt  number)  for  heating,  and 
a  very  small  and  possibly  negligible  increase  for  cooling. 

It  is,  of  course,  the  objective  of  the  present  work  to  resolve 
these  uncertainties,  since  it  is  felt  that  it  is  only  for  tur¬ 
bulent  heating  that  the  problem  is  well  in  hand. 

Both  the  work  of  Sze  and  that  of  Diessler  are  concerned 
with  the  variation  of  properties  in  the  radial  direction  only. 
The  effects  of  property  variation  in  the  axial  direction  have 
apparently  not  been  investigated  analytically.  It  is  usually 
tacitly  assumed  that  property  variation  in  the  axial  direction 
has  little  or  no  effect  on  the  local  heat  transfer  character¬ 
istics.  Thus  the  local  surface  conductance,  for  example,  is 
typically  evaluated  using  the  local  transport  properties  and 
temperature  ratios.  One  of  the  objectives  of  this  work  is  to 
help  cast  some  light  on  the  validity  of  this  assumption. 
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o  -  CIRCULAR  TUBE 

□  -  FLOW  BETWEEN  PARALLEL  PLANE3 


SECTION  IV 

EXPERIMENTAL  APPARATUS 


Tne  experimental  apparatus  used  in  this  work  consists 
of  three  elements:  l)  an  air  flow  metering  and  control  sys¬ 
tem;  2)  a  test  section  for  heating  air;  and  3)  a  test  section 
for  cooling  the  air.  These  elements  are  shown  in  Pigs.  2,  3, 
4,  and  5. 

The  air  flow  metering  system  is  shown  in  Fig.  1.  The 
air  flow  is  maintained  constant  by  the  use  of  a  series  of 
pressure  regulators  and  valves.  The  compressor  begins  oper¬ 
ating  when  the  tank  pressure  falls  to  80  psi  and  cuts  out 
when  the  pressure  reaches  100  psi.  A  regulator  on  the  outlet 
of  the  storage  tank  further  reduces  the  pressure  to  40  psi, 
and  another  regulator  in  series  again  reduces  the  pressure  to 
about  10  psi.  As  the  flow  rates  involved  in  the  laminar  flow 
work  are  small,  it  is  necessary  to  bleed  off  seme  air  between 
the  last  pressure  regulator  and  the  flow  meleis.  In  this  way 
it  is  possible  to  completely  avoid  any  detectable  variation 
in  the  air  flow  rate. 

Two  Fischer  and  Porter  Rotameters  are  used  to  measure 
the  flow  rate,  one  for  the  small  flows  needed  for  the  laminar 
work,  and  one  for  the  larger  flow  rates  needed  for  the  turbu¬ 
lent  work.  Both  were  carefully  calibrated,  and  the  experi¬ 
mental  uncertainty  is  estimated  at  ±  1  per  cent. 

The  apparatus  used  for  the  laminar  heating  experiments 
is  shown  in  Figs.  3  and  5.  It  consists  essentially  of  a  66 
inch  length  of  0.375"  0D,  0.005"  wall  nichrome  tube.  The 
first  30  inches  vre  unheated;  this  section  is  simply  a  de¬ 
veloping  section  for  the  velocity  profile.  The  next  30-inch 
section  is  heated  by  passing  an  electric  current  through  the 
tube  wall.  Thermocouples  are  spaced  along  tills  section  for 
the  measurement  of  the  tube  wall  temperature,  Fig.  5.  Be¬ 
cause  of  the  thinness  of  the  nichrome  tube  vail,  stem  con¬ 
duction  in  the  thermocouples  could  have  been  a  serious  prob¬ 
lem.  To  prevent  this  as  much  as  possible,  the  thermocouple 
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leads  are  not  taken  directly  away  i'rom  the  tube,  but  rather 
are  first  wound  around  the  tube.  The  heated  section  was  in¬ 
sulated  with  Johns-Manville  high  temperature  Cil-O-Cel  and 
magnesia  sheet  insulation. 

The  point  at  which  heat  transfer  begins  is  localized  as 
much  as  possible  by  introducing  the  electric  currert  through 
a  thin  flange  welded  to  the  tube,  as  shown  in  the  figures. 

It  is  not  necessary  to  know  precisely  where  the  heating  sec¬ 
tion  ends,  and  because  of  this,  and  because  of  the  high  tem¬ 
peratures  encountered  at  the  downstream  end  of  the  heating 
section,  a  simple  nichrome  collar  was  fitted  to  the  tube  for 
the  attachment  of  the  power  leads. 

Electric  power  was  supplied  to  *  he  nichrome  tube  from  a 
shop-built  power-pack  constructed  from  welding  transformers, 
with  a  variable  transformer  on  the  primary  for  control.  The 
current  to  the  tube  was  measured  using  a  Weston  current  trans¬ 
former  and  a  Weston  ammeter.  The  voltage  drop  along  the  tube 
was  measured  by  a  Ballantine  vacuum  tube  voltmeter  connected 
successively  to  a  series  of  voltage  taps  welded  along  the 
tube . 

For  purposes  of  making  energy  balances,  and  for  checking 
the  calculated  fluid  mean  temperature,  a  device  was  fitted  to 
the  exit  of  the  nichrome  tube  for  the  measurement  of  the  mean 
air  temperature  at  this  point,  see  Fig.  2.  A  ceramic  tube  is 
fitted  over  the  exit  end  of  the  nichrome  tube.  Short  lengths 
of  wire  traverse  the  ceramic  tube  and  serve  to  mix  the  air. 

At  the  end  of  this  ceramic  tube  a  1/8"  nozzle  is  fitted  to 
accelerate  the  flow  over  the  thermocouple.  To  reduce  the  heat 
losses  between  the  end  of  the  heating  section  and  the  exit 
temperature  thermocouple,  the  air  was  passed  back  past  the 
ceramic  tube  before  discharging  to  the  atmosphere. 

The  apparatus  for  the  cooling  experiments  is  somewhat 
more  complicated,  and  is  shown  in  Fig.  4.  From  the  flow- 
rater  the  air  flows  into  a  heater,  where  it  is  forced  through 
electrically  heated  nichrome  screens.  The  air  then  passes 
through  a  sharp-edged  orifice  where  it  is  turbulently  mixed. 
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A  thermocouple  is  situated  downstream  from  this  orifice  to 
measure  the  air  temperature.  The  air  then  passes  into  a 
.500"  OD,  .010"  wall  stainless  steel  tube.  The  first.  32 
inches  comprise  the  velocity  development  section.  Because 
of  the  low  air  flow  rates  involved  in  the  laminar  flow  ex¬ 
periments,  a  guard  heater  of  nichrome  ribbon  is  placed  around 
this  section  for  the  maintenance  of  the  air  temperature. 

Several  thermocouples  are  attached  to  trie  tube  wall  for  check¬ 
ing  this  temperature. 

From  tne  entrance  section  the  air  passes  through  a  con¬ 
stant  wall  temperature  section.  The  wall  temperature  is 
maintained  constant  by  passing  a  large  flow  of  water  over 
the  outer  surface  of  the  tube  in  a  counter-flow  heat  exchanger. 

The  exit  air  temperature  for  the  laminar  flow  experiments 
is  measured  with  the  thermocouple  system  shown  in  the  insert 
in  Fig.  <4.  A  1/8  inch  copper  plug  is  tightly  fitted  in  a 
balsa-wood  sleeve  which  in  turn  is  tightly  fitted  inside  the 
nichrome  tube  just  at  the  end  of  the  cooling  section.  The 
copper  plug  has  four  equi-spaced  No.  58  drill  size  holes  pass¬ 
ing  through  it  and  a  thermocouple  attached  at  its  downstream 
end. 

For  measuring  the  exit  air  temperature  during  the  tur¬ 
bulent  flow  experiments,  where  air  velocities  were  substan¬ 
tially  higher,  a  simpler  device  was  used.  This  consisted  of 
a  wire  mixer  and  a  l/8  inch  nozzle  to  accelerate  the  flow 
past  a  thermocouple. 

All  the  thermocouples  used  in  this  work  are  cromel- 
alumel.  The  thermocouples  are  connected  to  a  zone  box  which 
has  been  thermally  isolated  by  insulation,  and  in  which  cop¬ 
per  plates  hove  been  placed  to  insure  an  isothermal  zone. 

The  zone  box  connects  to  a  Leeds  and  Northrup  2-4-point  ther¬ 
mocouple  switch,  and  finally  to  a  Leeds  and  Northrup  Portable 
Precision  Potentiometer. 

A  summary  of  the  important  dimensions  of  both  the  heat¬ 
ing  and  cooling  apparatus  is  given  in  Table  I. 
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TABLE  I 

A  SUMMARY  OF  IMPORTANT  APPARATUS  DIMENSIONS 


Dimension 


Overall  tube  length 

Tube  inside  diameter 

Tube  wall  thickness 

Entry  length 

Entry  length  L/D 

Heat  transfer  section  length 

Heat  transfer  section  L/D 

Total  heat  transfer  area 


Gas 

Heating 

Appara  tus 

Gas 

Cooling 

Apparatus 

# 

66" 

62" 

C. 375" 

O.dBC" 

0.005" 

0.010" 

30" 

32" 

80 

66.7 

30" 

30" 

80 

62.5 

35.3  in£ 

45.3  in 

Includes  6  inches  extending  beyond  and  at  heating  section. 


18 


GA3  COOLING  APPARATUS 


SECTION  V 

REDUCTION  OP  EXPERIMENTAL  DATA  AND  TEST  RESULTS 

I'Hnlr.u?7  Heating  Experiments  --  In  the  laminar  flow  heat¬ 
ing  experiments  the  measured  quantities  were  the  tube-wall 
temperature  as  a  function  of  axial  position,  the  air  flow 
rate,  the  inlet  air  temperature,  and  the  total  power  input. 
With  this  information,  and  knowledge  of  the  magnitude  of  the 
heat  losses  through  the  surrounding  insulation,  the  local 
heat  transfer  rates  and  local  Nusselt  numbers  can  be  calcu- 
la  ted . 

Before  describing  the  actual  calculation  procedure,  the 
problem  of  heat  losses  through  the  insulation  should  be  dis¬ 
cussed.  With  the  very  low  air  flow  rates  involved,  on  the 
order  of  1 .  pounds  per  hour ,  this  problem  was  a  very  serious 
one.  Actually  the  magnitude  of  the  heat  transfer  through 
the  insulation  was  of  the  same  order  as  that  to  the  air  flow¬ 
ing  inside  the  ni chrome  tuie.  Increasing  the  amount  of  in¬ 
sulation  would  have  helped  very  little,  in  view  of  the  al¬ 
ready  small  tute-diameter-to-insulation-diameter  ratio.  The 
advantage  of  the  small  decrease  in  heat  losses  obtained  in 
this  way  would  have  teen  more  than  compensated  for  by  the 
increase  in  apparatus  response  time  which  would  accompany 
any  such  solution. 

Since  it  is  necessary,  in  order  to  calculate  the  local 
Nusselt  number,  to  have  accurate  knowledge  of  the  local  heat 
input,  it  was  necessary  to  know  the  magnitude  of  the  local 
heat  losses  accurately.  For  this  reason  a  number  of  test 
runs  were  performed  with  no  air  flowing.  From  the  measured 
axial  distribution  of  wall  temperature  and  the  total  power 
input,  it  was  possible  to  construct  a  chart  giving  the  local 
heat  loss  as  a  function  of  local  wall-temperature-to-ambient- 
temperature  difference  and  position. 

The  question  may  arise  as  to  whether  or  not  the  local 
heat  losses  obtained  in  this  way  are  the  same  as  those  which 
occur  during  the  actual  test  runs,  in  view  of  the  fact  that 
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the  axial  temperature  distributions  were  different.  A3  long 
as  the  major  share  of  the  thermal  resistance  is  in  the  insu¬ 
lation,  and  as  long  as  the  first  derivative  of  the  temperature 
versus  axial  position  curve  is  small,  the  local  heat  losses 
as  obtained  previously  should  correspond  to  the  losses  under 
the  actual  test  conditions.  These  conditions  are  least  well 
satisfied  in  the  entry  to  and  the  exit  from  the  test  section. 
Consequently  the  results  from  these  regions  contain  greater 
uncertainty  chan  do  the  results  from  other  regions  of  trie 
test  section.  Because  of  the  possible  error  in  the  exit  re¬ 
gion,  the  calculations  were  not  carried  into  this  region. 
Calcula cions  for  the  great  bulk  of  the  test  section  exclud¬ 
ing  this  exit  region,  using  the  axial  test  run  temperature 
distributions,  indicate  that  the  conditions  enumerated  above 
are  adequately  satisfied  and  that  the  use  of  the  heat  loss 

e* 

data  obtained  under  conditions  of  no  air  flow  introduces 
negligible  error  in  the  final  result. 

With  heat  loss  as  a  function  of  axial  position  and  tem¬ 
perature  determined,  the  calculation  procedure  is  as  follows. 
Using  the  measured  wall  temperature  distribution,  the  local 
heat  loss  as  a  function  of  position  is  determined. 

L  =  L(x)  (16) 


From  the  power  input  measurements  and  the  local  heat  loss 
the  local  rate  of  heat  transfer  to  the  air  flowing  inside 
the  tube  can  be  found.  Thus, 

q'  =  P(x)  -  L(x)  =  q ' (x)  (17 ) 

Check  calculations  indicate  that  for  Runs  1  to  4  inclusive 
the  local  power  input  can  be  assumed  to  be  the  same  as  the 
average  power  input  over  the  length  of  the  test  section. 
Excluding  the  first  2  inches  of  the  test  section  where  the 
unct  tainty  due  to  the  uncertainty  in  L(x)  is  high,  the 
effect  of  this  assumption  is  to  change  the  calculated 
resul;  by  less  than  1  per  cent.  For  these  runs 

P  ( x )  =  P 
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Having  now  the  heat  transfer  rate  to  the  air  inside  the 
tube  os  a  function  of  position,  the  local  mixed  mean  air  tem¬ 
perature  can  be  determined.  This  is  obtained  by  a  simple 
energy  balance  on  the  volume  Included  within  the  tube  from 
the  start  of  the  heating  section  to  the  local  point  of  in¬ 
terest.  The  resulting  equation  is 


/  q'dx 
0 


(16) 


1 


It  is  now  possible  to  calculate  the  local  surface  con¬ 
ductance.  From  the  defining  equation  for  the  surface  con¬ 
ductance 

h  =  q'7(Tw  -  Tm)  (19) 

The  rate  of  heat  transfer  per  unit  area  is  given  by 


q"  =  q'/27rro 


Substituting  Eqs.  (2C)  and  (l8)  into  (19)  gives 


h(x) 


P  -  Mx) 

2irr 

o 


X 

/  q'dx 

0  J 


(20) 


(21) 


The  local  Nusselt  number'  can  then  be  determined  from 

NNu  =  bD/k  (22) 

Actually  the  calculation  procedure  was  somewhat  more 
complex  due  to  the  fact  that  the  air  properties  which  en¬ 
tered  into  the  above  calculations  for  h(x)  are  themselves 
functions  of  temperature.  In  effect,  a  local  mean  tempera¬ 
ture  was  first  assumed  and  used  to  evaluate  a  mean  c  for 

P 

the  range  from  TQ  to  Tm  .  Using  this  value  of  cD  a  new 

mean  temperature  could  be  calculated.  The  mean  c  from 

P 

T,j  to  Tm  was  then  evaluated  using  the  calculated  value  of 
T  ,  This  iteration  procedure  was  followed  until  no  signi¬ 
ficant  change  resulted. 

All  air  properties  used  in  these  calculations  were  taken 


23 


: 4. -  __1.‘ 'CVS 


erw  .friniiE; witwaj.  c*\i 


* 


I 


from  the  National  Bureau  of  Standards  Bulletin  No.  56^. ‘‘“k 
In  runs  5  and  6  the  variation  in  tube  electrical  re¬ 
sistance  makes  the  assumption  of  constant  power  Input  un¬ 
tenable,  and  the  NNu  calculated  using  thl3  assumption 
must  be  corrected  as  follows: 

M  _  hD  _  q1 
‘  Nu  k  kttAT 


Since  k  is  not  a  function  of  the  electrical  resistance, 

R  , 


dN 


Nu 


N 


Nu 


do 


q‘ 


dAT 


AT 


If,  as  is  the  case  nere,  Aq  '  «  q'  and  AAT  <<  AT 


AN 


Nu 


N 


Nu 


Aq '  AAT 
q '  AT 


where  Aq 1  and  AAT  are  the  changes  in  q'  and  AT  due 
to  the  variation  of  tube  resistance,  and  where  q1  and  AT 
are  the  values  computed  previously.  Then, 


P ( x )  =  IdR(x) 


Thus 


AP'  =  r-AR(x) 
where  Aft(x)  is  given  by 

AR(x)  =  Rabt 

where  a  is  the  temperature  coefficient  of  resistance  of 
the  tube,  R  is  the  length  mean  resistance,  and  where 

6t  =  Vx>  '  T.  /  T„(x)dx 


0 


W 


Thus  the  actual  power  input  is  given  by 

P(x)  =  i2r  +  i2ar 

=  i2r(i  +  ^|) 

P(x)  =  p(x)  [  1  +  Qt6t  ] 
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where  a  bar  over  a  quantity  is  used  to  indicate  its  cor¬ 
rected  value. 

Thu  3 

q'  =  P(x)  -  L(x) 

q'  =  P ( x ) [ 1  +  a6t]  -  L(x) 


Aq '  Is  easily  found  by  subtracting  q'  from  the  value  of 
q1  calculated  previously. 

AAT  can  be  found  as  follows: 

—  i  x  _ 

AT  =  T  -  T - —  /  a'dx 

u  wc  n 

p  0 
or 


Thus 


^ -  2  q'AL 

P 


AAT  =  AT  -  AT 


T  =  -Z  +  2  =  _2_  [2  qlAL  ..  2  q'AL] 

wc  wc  wc  1  1  M  J 


The  Nusselt  numbers  obtained  above  were  plotted  against 
a  nondimenslonal  distance  x+  ,  defined  earlier.  As  x+ 
involves  the  properties  of  the  air,  the  question  again  arises 
as  to  where  these  proper';!  es  should  be  evaluated.  It  was 
found  that  use  of  the  local  mixed-mean  temperature  brought 
about  the  test  agreement  between  the  experimental  arid  the 
theoretical  results. 

As  is  well  known,  the  analytical  solutions  for  the  local 
Nusselt  number  are  dependent  on  the  distribution  of  heat  in¬ 
put.  Since  the  local  heat  input  in  all  the  runs  was  very 
closely  approximated  by  a  step  plus  a  ramp  (negative),  that 
Is 

q"  a  qQ[l  +  bL]  (2j) 

the  results  of  a  calculation  for  this  type  of  heat  Input 
should  apply  to  all  the  runs.  It  should  be  noted,  however, 
that  In  general  the  value  of  b  was  different  for  different 
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run3.  Fro:.  general  solution  for  variable  heat  Inpu 
(see  Section  II) 

r  x+ 

T  -  T  =  rr  !  K (x  -  L,:'+)q(L)dL 

o  ,<  o 

or  for  rf  =  1 


w  o  k 

Substituting  for  g(x+  -  L,l) 

+ 


T--  -  T.  =  /  g(x+  -  L,l)q(L)dL 


0 


r  x 

T  -  T  =  —  f 

w  0  '•  n 

K  0 


4  - 


-y2(x+-L) 

e  m'  ' 

S  ~T~~,  5T” 


q(L  )aL 


For  a  step-plus-a -ramp  type  heat  input, 

q(L) 


q”(l  +  tL) 


Substituting, 


r  x 
T  -  T  =  q"  —  f 

T.r  r\  ^  r\  J 


w  o 


k  0 


4  - 


z 

rr. 


e"''m<x+-L> 


[1  I  tLjdL 


Integrating  between  the  given  limits. 


'w  o  -  Mo 


k 


f. 

x+  cx+" 

t 

A  i- 

2  J 

+  Z 


2  + 
e"Kmx 

•/V  (-y2) 

rm  '  rm' 


-bZ 


2  + 
e-W 

m  y2h'(-y2) 
rm  v  'm; 


2  + 

e^m*  ( 

4+ 

2 

-v) 

' 

- 

J 

Also,  from  energy  considerations. 


.+ 


4r  x" 

Tm  -  T  =  —2.  f  q"(x+)dx+ 
m  0  k  C 


-  Z 

rr. 


”ZJ - 

Y  H 1 

■  m 


m 


Tm  *  To  -  <  l  U  +  bx+>dx+ 


4r 


T  -  T  = 
m  o 


+  bx 
x  +  — 


+2 


£0 


Combining  yields 


T  -  T  a  a"  —  <  2-7 

w  -  *o  .  1  rn  J 


-v  x 
rrr. 


m  7  H' (-v  )  n  V  H' 1  -y  ) 
'm  '  'nr  ^m  k  'nr 


-y  x 
'  m 


:r- 

<  •»  in 


x  1  W  1 
—  -  IT  r  ~T 
y  y  y 

1  m  '  m  r  in 


Since  Nn  =  hD/k 


x  w  rn 


2r 

o 

nr 


(T  -  T  )  k 
q”(l  +  bx+)  r{ 


2  + 
-Y  x 
fm 


2  '  1  +  tx  m  y>(-y2) 


£  2„,  ,  27  2  ZT  + 


2  + 
-Ymx 


Simplifying  further 


y  W  (-7  ) 

rm  v  fm' 


-  32  g  rr 


From  the  arbitrary  heat  input  solution. 


G  11 

52  2  -£  =  — 

n  A  24 
n 


Therefore , 


1  ]  1 

-2  — r  ~  = 

m  y  H'(-y^)  24 

'  m  '  ■  m  ' 


*■  «-  X. 


Therefore, 


1  1  1 


-yV 

rmA 


X  2 1  +  bx+  m  ^777  + 7 


y>(<)  7 


-?y 


Define 


A(x  +  )  = 


'7X  + 


&'(•&  +  24 


B(x+)  =  Z 


rr.  .  ,.c 


Tnen 


r„H'(-y  )  v 2 

m  '  rrri'  'rr. 


X-  e'^  + 


2v  \~2  + 


Nux  r^7  [A(x d  -  tB^+) 

or 

2(1  +  fcx+) 

wNu  =  — : j: - - — — 

x  A(x  )  -  bB(x+) 

Values  of  A(x  +  )  and  B(x  +  )  ha7P  v 

r.A .  _  .  'x  ‘  have  teen  evaiua 


(2-f ) 


(25) 


(26) 


given  in  Table  II. 


n  evaluated  and  are 


TABLE  II 


THE  FUNCTIONS  A(x*)  AND  -B(x+) 


+ 

X 

A(x+) 

-B(x+)  x  10*5 

0.02 

.5232 

4.919 

0.04 

.3843 

12.09 

o 

• 

o 

cm 

.4156 

20.14 

0.08 

.4  329 

28.61 

0.10 

.4431 

37.36 

0.12 

.4492 

46.52 

0.14 

.^529 

56.26 

0.16 

.4551 

64.55 

The  values  of 

b  used  In  Eq. 

(26)  are  found  from  the 

heat  Input  given  by  Eq.  (23). 


Cooling  Experiments  —  In  the  cooling  experiments,  the 
local  Nusselt  numbers  were  not  determined;  only  the  length- 
mean  Nusselt  numbers  were  calculated.  The  experimental  data 
taken  consisted  of  the  air  flow  rate,  several  temperatures 
along  the  entry  length,  the  water  inlet  and  outlet  tempera¬ 
ture,  the  air  exit  temperature,  and  the  water  flow  rate. 

Again  the  evaluation  of  both  the  experimental  value  of 
the  Nusselt  number  and  the  theoretical  solutions  Involve  the 
problem  of  what  temperature  to  use  when  evaluating  the  prop¬ 
erties.  Since  the  data  Is  for  the  length-mean  Nusselt  num¬ 
ber,  it  Is  necessary  to  use  some  sort  of  length-mean  tem¬ 
perature.  If  the  thermal  transport  properties  are  constant, 
and  If  the  surface  conductance  (h)  Is  constant,  the  length- 
mean  temperature  is  the  simple  log-mean  temperature  which  Is 
usually  employed  In  heat  exchanger  work.  In  the  actual  sit¬ 
uation,  however,  this  Is  rot  a  true  length-mean  temperature 
for  two  reasons.  First,  the  thermal  conductivity  Is  not 
constant;  actually  k  is  larger  at  the  front  part  of  the 


cooling’  section  where  the  air  is  at  a  higher  temperature; 
and  secondly,  since  the  thermal  boundary  layer  is  not  fully 
developed  at  the*  entrance  to  the  cooling  sect! or  ,  the  Mussel t 
number  there  is  considerably  higher  than  at  the  exit. 

In  the  case  of  the  turbulent  data  the  mean  temperature 
was  taken  'o  be  the  simple  log-mean  temperature,  as  this  is 
by  far  the  easiest  and  since  the  entry  lengtn  effect  is 
small . 


In  the  laminar  case  the  uuta  was  reduce:  tnree  times; 
once  with  the  properties  evaluated  at  t re  simple  log-mear 
temperature,  once  with  the  properties  evaluated  at  the  exit 
temperature,  ar.u  once  with  the  properties  evaluated  at  a 
mean  temperature  calculated  from  a  Graetz  solution  which  as¬ 
sumes  constant  properties  but  not  fully  developed  thermal 
boundary  layer.  Tine  evaluation  at  the  log-mean  temperature 
was  done  primarily  for  general  interest  and  also  to  demon¬ 
strate  the  sensitivity  of  the  results  to  the  choice  of  tem¬ 
perature  at  which  the  properties  are  evaluated.  The  actual 
length-mean  temperature  will  lie  between  the  length-mean 
temperature  given  by  the  Graetz  solution  and  the  exit  tem¬ 
perature,  for  the  reasons  discussed  previously.  Hence  the 
effects  of  the  radial  variation  of  properties  should  lie  be¬ 
tween  the  effect  indicated  when  the  properties  are  evaluated 
at  the  Graetz  solution  mean  temperature  and  the  effect  indi¬ 
cated  when  the  properties  are  evaluated  at  the  exit  temperature. 

The  actual  data  reduction  procedure  for  the  turbulent 
data  was  as  follows:  Since  the  water  temperature  was  con¬ 
stant  within  a  few  degrees,  the  log-mean  temperature  can  be 
calculated  from 


T 

w 


+ 


AT.  -  AT  . 
_ in _ out 

in  AT,  /AT  , 
in'  out 


(27) 


The  effectiveness  can  be  evaluated  from 


e 


T 


h 


In 


m 

^OUt 


T^in  "  Tcout 


(28) 
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7^1  F  —  l  -  )°-6(2r  )  (53) 
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Except  in  the  initial  runs  where  the  water,  flow  rate  was 
kept  low  so  as  to  have  a  significant  water  temperature,  change 
with  which  to  make  an  energy  balance  for  check  purposes,  the 
thermal  resistance  of  the  water  side  was  negligible,  i.e., 
less  than  1  per  cent  of  the  total. 

The  length-mean  Nusselt  number  can  then  be  evaluated 

from 

NNu  =  hD/k  (22> 

where  k  is  evaluated  at  the  log-mean  temperature  calcu¬ 
lated  above.  The  Reynolds  number  can  be  obtained  from 
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He 


where  p.  In  evaluated  at  the  log-mean  temperature .  T:.e  re¬ 
sults  we  re  then  norm  11  zed  by  dividing  tie  experiment  a  1 
Mussel  t  number  by  're  Mussel*,  rurr.fer  giver  t  y  '  :.e  soldi!  e* 
of  Sleicr.er  a:  J  Trirus  using  the  N.  ottair.cd  above  tn 

ac 

evaluate  trie  exi*  x+  . 


For  sons'. a:  t  wall 
given  by 
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‘Nu„  , 
•iso 


Gx  + 


G  .  1  +• 
87.  Ui  e'^r‘x 


This 


Nu 


has  been  evaluated  for  various  values  of 


For  the  apparatus  in  use  x/rQ  Is  fixes.  Assuminr 
N„  =  constant, 


Pr 


iso 


"  ^  Re  ^ 


’He  * 


Ir  tr.e  lir-.ir.ar  case,  using  the  log-mean  temperatures, 

the  data  reduction  procedure  was  the  same,  except  of  course, 

that  the  appropriate  eigenvalues  ana  constants  from  Appendix 

A  were  used  far  the  calculatior.  of  Nm,,  . 

^ui  so 

For  the  cecond  method  of  reduction,  where  the  properties 
’"ere  evaluate;  at  the  mean  temperature  given  by  a  Graetz  so¬ 
lution,  the  calculation  procedure  was  somewhat  more  complex. 
First  the  log-mean  temperature  was  evaluated  as  above,  and 

4- 

also  an  exit  x  .  From  the  Graetz  solution  a  nondi mensional 

mean  temperature  can  be  obtained  as  a  function  of  the  exit 

x+  (see  below).  Tnls  new  mean  temperature  is  used  to  evalu- 

+ 

ate  a  new  exit  x  ,  and  the  procedure  repeated  until  the 
value  of  the  mean  temperature  is  stable.  A  Nusselt  number 
based  on  the  properties  at  this  temperature  is  then  obtained. 
This  Nusselt  number  is  then  normalized  by  divining  It  by  the 
mean  Nusselt  number  Njju,  predicted  by  the  Graetz  solution 
for  the  final  value  of  x+  . 

A  mean  temperature  based  on  the  Graetz  solution  can  be 
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obtained  as  follows:  From  Eq.  (y)  the  local  mixed-mean 
temperature  at  any  xf  is  give:  by 

»  2  + 
jr  -X,  x 

3z  e  “ 
n  K 

The  length-mean  temperature  is  given  by 


e„,  = 

til 


0  1 X 


+ 


Substituting  an: 


r.tegra  t 


;£ 


yiei 


c- 

o 


•  xt. 


■4“ 


ane 
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Analysis  of  Exper imenta  1  Uncertainty  —  An  analysis  of 
tre  experimental  uncertainty  is  presented  in  Appendix  C.  It 
is  concludes  that,  employing  the  above  described  methods  of 
reduction  of  data,  the  uncertainty  in  the  measured  Nusselt 
numbers  for  tre  laminar  .seating  experiments  is  approximately 
±  11  per  cent.  For  the  laminar  coolirg  experiments  an  un¬ 
certainty  of  +  p  per  cent  in  Nusselt  number  Is  estimated, 
and  for  tne  turbulent  cooling  experiments  the  uncertainty 
in  Nusselt  number  is  estimated  at  +  6T  per  cent. 

Test  Results  --  The  complete  experimental  data  for  lami¬ 
nar  heating,  laminar  cooling,  and  turbulent  cooling  are  pre¬ 
sented  in  tabular  form  in  Appendix  B. 

The  laminar  heating  experiments  consist  of  a  series  of 
six  test  runs  each  with  an  entering  Reynolds  number  near 
1200.  Data  were  then  recorded  at  a  number  of  positions  along 
the  tube.  It  will  be  noted  that  the  local  Reynolds  number 
decreases  along  the  tube  because  of  the  increase  nf  viscos¬ 
ity,  this  effect  being  most  marked  for  the  high  temperature 
runs.  Each  run  represents  a  successively  higher  heating 
rate  with  wall  surface  temperature  reaching  1100°F  for  Run  6, 
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and  the  nn-'face-to-fluid  temperature  difference  reaching 


-*C  n 

*;>  ( 0  r  • 


test  runs  a'  varying  Reynolds  numbers  wit:,  the  air  inlet 
ter  j  eraturc*  successively  increase::  to  a  tr.axlr.urr.  of  37 *3°F . 

T:.e  turculer.t  cooling  experimer Is  cover  a  Reynolds 
nur.ter  range  fror.  about  1)000  to  )5000 ,  with  air  inlet  tem¬ 
perature  varylrg  up  to  a  maxi  nun.  of  1N00°F. 

T.ne  test  results  are  presented  graphically  on  Kirs,  o 
to  1  Tr.e  laminar  heating  data  are  shown  on  Figs .  6  to  11 
where  the  local  Nusselt  .number  is  plotted  as  a  furctio:  of 
the  nondimen.sional  lengti  parameter,  x+  .  On  Fig.  12  t .-  e 
ratio  of  local  measured  Nusselt  number  to  predicted  Nusselt 
number  is  plotted  as  a  function  of  the  absolute  wall-to-mcar- 
fluid  temperature  ratio.  The  laminar  cooling  data  are  simi¬ 
larly  plotted  or.  Figs.  1),  Id,  and  15,  using  t-.e  three  meth¬ 
ods  of  evaluating  fluid  properties  described  previously.  The 
turbulent  cooling  data  are  shown  or.  Figs.  16  and  17,  Or1 
Fig.  16  the  measured  overall  mean  Nusselt  numbers  are  plotted 
as  a  function  of  Reynolds  number.  On  Fig.  17  the  ratio  of 
Nusselt  rumrer  to  the  Nusselt  number  predicted  by  tr.e  Sleicner 
and  Tribus  solution  is  plotted  as  a  function  of  the  tempera¬ 
ture  ratio. 


i  -  .u.  -siry 
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SEC 710V  VI 

DISCUSSION'  OF  TEST  RESULTS 


Laminar  Heating  —  The  experimental  values  of  the  local 
Nvu  plotted  against  the  local  values  of  x+  are  shown  in 
Figs.  6  to  11  inclusive.  Also  plotted  as  a  dashed  line  are 


the  values  of  tne  local 


W 


Nu 


as  predicted  by  the  constant 


It  should  appear  on  a  plot  of  ^Nu'^Nu-' s  vs 


T  /T 
w  m 


Nn  are  evaluated  at 


property  analytical  solutions  for  the  actual  heat  rates  which 
existed  In  each  case. 

It  Is  seen  that  the  experimental  values  correspond  well 
to  the  tneoret ic.- 1  solutions,  even  when  the  temperature  dif¬ 
ferences  involved  are  large  and  the  thermal  profile  is  de¬ 
veloping.  Tne  agreement  is  actually  test  in  the  runs  where 
the  temperature  differences  involved  were  highest. 

If  there  were  an  appreciable  variable  properties  effect 

,  Fig. 

12.  The  Lest  straight  line  correlation  indicates  that  the 
variable  properties  effect  is  negligible  if  the  properties 
used  in  tne  evaluation  of  the  local 
the  local  mixed-mean  temperature.  Because  of  the  scatter  in 
the  experimental  data,  and  the  rather  large  estimated  experi¬ 
mental  uncertainty,  a  small  variable  properties  effect  could 
be  masked,  rut  fh^re  is  certainly  no  substantial  effect,  at 
least  for  temperature  ratios  up  to  2.0  .  This  is  at  vari¬ 
ance  with  the  analytical  solutions  which  predict  a  drop  in 
conductance  for  heating.  However,  the  analytical  solutions 
do  not  take  Into  consideration  any  effects  of  the  develop¬ 
ment  of  the  temperature  profile,  or  the  effects  of  property 
variation  with  tube  length,  3rd  it.  Tull y  well  be  that  this  is 
an  effect  that  compensates  for  the  Influence  of  the  distor¬ 
tion  of  the  velocity  and  temperature  profile  at  any  particu¬ 
lar  cross-section. 

Laminar  Cooling  --  The  experimental  values  of  N,T  /it, 
- a  r  Nuj_so 

vs  Tm/Tw  for  the  laminar  flow  cooling  experiments  are 
shown  in  Figs.  1 3,  14,  and  id.  In  Fig.  13  all  properties 
were  evaluated  at  the  log-mean  temperature.  In  Fig.  14  all 
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properties  were  evaluated  at  the  mean  temperature  given  by  a 
Graets  solution,  as  discussed  earlier.  In  Fig.  1*3  all  prop¬ 
erties  were  evaluated  at  the  exi*  temperature. 

In  all  these  presentations  the  best  straight  line  cor¬ 
relation  intersects  the  Nu/Nu,  axis  at  about  0.9;<,  in¬ 
dicating  tnat  the  experimental  results  are  ahcut  six  per  cent 
low.  This  is  consistent  with  the  estimated  experimental  un¬ 
certainty,  and  probably  means  that  the  major  sources  of  un¬ 
certainty  are  fixed  rather  than  random. 

As  was  discussed  in  the  section  or,  data  reduction,  rone 
of  the  temperatures  used  for  the  evaluation  of  the  proper¬ 
ties  ore  the  actual  length-mean  temperature.  It  is  to  he  ex- 
pectec.  that  the  true  length-mean  temperature  lies  somewhere 
between  the  "Graetz  mean  temperature"  and  the  exit  tempera¬ 
ture.  Hence  the  true  value  of  "a"  should  be  somewhere  be¬ 
tween  the  value  of  the  slope  of  the  straight  line  correlation 
in  Figs.  14  and  15.  It  is  seen  that  Sze 1 s  predicted  value 
of  a  =  -0.0&  does  lie  within  this  range.  Hence,  within  the 
experimental  uncertainty  it  would  appear  Dial  Sze 1 s  predic¬ 
tion  for  the  case  of  laminar  cooling  does  apply.  However,  it 
should  be  noted  that  this  Is  a  quite  small  correction,  and 
the  important  conclusion  that  one  can  draw  is  that  both  an¬ 
alysis  and  experiment  are  consistent  in  showing  a  small  ef¬ 
fect  io  tiie  same  airection. 

In  Fig.  13,  v/ here  the  properties  were  evaluated  at  the 
log-mean  '"emperature ,  the  temperature  dependent  properties 
effect  appears  to  be  negligible  (a  =0).  Thus  it  would  ap¬ 
pear  that  if  the  properties  are  evaluated  at  the  log-mean 
temperature  the  variable  properties  effect  could  be  ignored, 
in  that  it  would  be  accounted  for  by  evaluating  the  proper¬ 
ties  at  the  log-mean  temperature.  This  Is  not  entirely  cor¬ 
rect  for  the  following  reason.  With  the  particular  apparatus 
used  for  this  work,  trie  effects  of  property  variation  were 
just  compensated  by  the  error  in  evaluating  the  properties 
at  the  log-mean  temperature.  In  general,  however,  if  the 
geometry  of  the  apparatus  were  changed  one  would  not  expect 
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this  same  result  because  It  would  then  be  possible  to  dupli¬ 
cate  the  same  values  of  ”„/T.  (T„  «  log-mean  temoera lure ) 

while  altering  the  actual  length-mean  temperature  and  thus 
the  variable  properties  effect.  Again,  nowever,  since  tne 
effect  is  evidently  small,  it  would  seem  tnat  a  satisfactory 


design  procedure  would  ce  to  evaluate  all  properties  at  log- 
mean  temperature  a:.:  Ignore  any  effect.'  of  temperature  de- 
pender.t  properties. 


Tsrtulent  Cooling  —  Tne  experimental  values  of  NMt, 

w  *  NU 

are  shown  as  a  function  of  NRe  in  Fig.  16.  It  is  seen  tha 
the  data  points  are  consistently  above  the  values  predicted 
by  the  commonly  used  lully  developed  turbulent  flow  correla¬ 
tions.  This  is  partially  due  to  tne  entrance  effect;  how¬ 
ever,  even  if  this  is  accounted  for,  the  data  points  are 
stili  cot oister.tly  about  10  to  20  per  cent  higt.  Experimen¬ 
tal  error  might  be  suspected,  if  all  other  data  was  consist¬ 
ent  witr.  the  common  correlations.  However ,  the  experimental 
values  of  i\’^u  cor.  csgui.i  very  .'.ell  to  tne  Sleicher  and 
Tribus-  solution.  Their  solution  in  turn  was  based  on  their 
own  experimental  ’work.  It  Is  not  known  if  any  explanation 
for  the  discrepancy  between  the  results  of  Sleicher  and 
Tribus  and  tne  common  correlations  exists. 

The  experimental  values  of  NNu/NNul£^  (where  NNu1sq 
is  the  value  predicted  by  Sleicher  and  Tribus)  plotted 
against  Tm/Tw  are  shown  in  Fig.  17.  As  can  be  seen,  the 
data  is  best  represented  by  “  1.00  .  Despite 

the  discrepancy  between  the  present  data  and  the  common  cor¬ 
relations  (or  the  discrepancy  between  the  Sleicher  and  Tribus 
solution  and  the  common  correlations),  it  seems  evident  that 
there  is  negligible  temperature  dependent  properties  effect 
over  the  temperature  ratio  range  1.5  to  3.0,  provided  that 
all  fluid  properties  are  evaluated  at  tire  meat  temperature 
will  'especl  to  length  (the  log-mean  in  this  case). 
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SECTION  VII 


CONCLUSIONS 

The  conclusions  of  this  I nvest Iga tior  -nay  be  summarize  ! 
as  follows: 

1.  Two  conclusions  car.  he  reached  from  the  experimental 
data  for  laminar  heating.  The  first  is  t ha ’  the  Graetz  solu¬ 
tion  and  its  extensions  seem  to  preuict  correc‘ly  the  varia¬ 
tion  of  local  Nusselt  number  if  the  air  properties  are  evalu¬ 
ated  at  tne  local  rr.ixed-mean  temperature.  That  is  ”o  say, 
tne  thermal  entry  effect  seems  to  correspond  generally  to 
that  predicted  by  the  proper  analytical  solution,  ever,  though 
the  assumption  of  constant  properties  made  in  these  solutions 
does  not  obtain.  The  second  conclusion  is  that  the  actual 
effect  of  the  radial  variation  o'"  properties  on  the  local 
Nusselt  number  is  small,  at  least  for  temperature  ratios  up 
to  2.0  .  The  best  fitting  straight  line  correlation  (see 
Fig.  12)  indicates  that  this  effect  is  negligible  provided 
all  properties  are  evaluated  at  the  local  mixed-mean  tempera¬ 
ture.  Because  of  the  scatter  in  the  experimental  data  ar. i 
the  rather  large  experimental  uncertainty,  a  small  property 
variation  '■'O'l^d  b°  ^'idde,n,  but  cr>T1tair. ly  ar.  effect  as 

large  as  predicted  does  not  seem  to  obtain. 

2.  For  tne  cooling  of  a  gas  in  laminar' flow  the  da "a 

(j 

obtained  were  only  for  the  mean  Nusselt  number  with  respect 
to  tube  length,  rather  than  the  local  Nusselt  number,  but 
again  the  Graetz  solution  is  apparently  accurately  applicable 
if  all  fluid  properties  are  evaluated  at  a  correct  mean  fluid 
temperature  with  respect  to  tube  length.  For  constant  wall 
temperature  cooling,  the  logarithmic  mean  temperature  is  suf¬ 
ficiently  close  to  the  true  mean,  and  in  some  cases  is  actu¬ 
ally  a  better  temperature  at  which  to  evaluate  properties 
than  the  true  mean.  This  conclusion  is  based  on  absolute- 
mean-temperature-to-wall-temperature  ratios  up  to  2.0  . 

3.  For  a  gas  being  cooled  in  turbulent  flow,  no 


■vA- 1 I  j  <r  1  Vw.JV.l.aviivi 
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temperature  dependent  properties  effects  could 


sign! ficant 

te  detected  for  abaolule-rr.eun-temperature-to-wall-texperature 
ratios  from  1.5  to  3.0  .  The  experiments  covered  only  mean 
!’  "selt  numbers  with  respect  ‘o  tube  length  with  constant 
wall  temperature  cooling,  and  all  fluid  pru_per'ties  were  evalu¬ 
ated  at  the  logarithmic  mear  ,  wrier,  for  turbulent  flow  Is 
very  close  to  the  actual  mean  temperature  with  respect  to 
length. 

Tr is  cone  1 us lor  is  in  agreement  with  the  experimental 

,  "> 

results  of  Zell  nix  and  Churchill who  also  found  that  the 
usual  correlations  sold  for  tne  of  turbulent  cooling  if 

all  properties  are  evaluated  a"  the  mixed -me a:  temperature. 


u.  Consideration  of  the  four  significant  regimes  of 
neat  transfer  in  a  circular  tube,  laminar  ana  turbulent  : low, 
heating  ana  cooling,  leads  to  the  conclusion  t.nat  it  is  only 
in  tur t  alert  f low  heating  that  there  is  a  significant  tem¬ 
perature  dependent  properties  effect  for  gases,  a*  least  up 
to  absolute  temperature  ratios  of  2.0  or  ;.0  .  In  all  other 
cases  the  constant  property  analytical  solutio:  s,  or  tr.e 
small  temperature  difference  experimental  correlations,  are 
believed  to  be  of  acceptable  accuracy  for  engineering  calcu¬ 
lations  If  all  fluid  properties  are  evaluated  at  the  local 
mixed-mean  temperature  where  local  Nusselt  numbers  are  con¬ 
cerned,  and  the  mean  temperature  with  respect  to  length  in 
tne  case  of  mean  Nusselt  numbers.  For  turbulent  flow  heating 
the  same  procedure  should  be  follo’wed,  with  the  exception 
that  the  resulting  Nusselt  number  should  be  multiplied  by 
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APPENDIX  A 


SUMMARY  OF  EIGENVALUES  AMD  CONSTANTS 
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SUMMARY  OF  EXPERIMENTAL  RESULTS 
LAMINAR  HEATING  EXPERIMENTS 
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APPENDIX  C 

ANALYSIS  OF  EXPERIMENTAL  UNCERTAINTY 


Laminar  Heating  Experiments  --  The  uncertainty  in  the 
experimental  determination  of  Nusselt  number  for  the  laminar 
heating  experiments  can  be  estimated  by  an  analysis  of  the 
uncertainty  of  the  various  items  of  information  going  into 
Eqs.  (21)  and  (22).  It  will  be  noted  that  the  tube  diameter 
cancels  out  and  the  pertinent  variables  are  power  measure¬ 
ment,  heat  leak,  flow  rate,  the  temperature  difference  be¬ 
tween  the  tube  wa 1 1  and  the  entering  fluid,  and  the  fluid 
properties.  Two  test  runs,  representing  extremes  of  the 
test  data,  have  been  analyzed,  and  the  following  tables  list 
the  magnitude  and  estimated  probable  uncertainty  in  the  per¬ 
tinent  variables  (neglecting  the  specific  heat  for  which  the 
probable  uncertainty  is  small  relative  to  the  others). 

nun  No.  1,  Position  20" 

Variable  Magnitude  Uncertainty  Relative 

Interval  Uncertainty 

Gross  Power, 


Btu/in 

0.846 

— 

± 

0 . 01 

Heat  Leak 

Btu/ In 

0.460 

— 

+ 

C.04 

(Tw  -  V'  °F 

50.5 

±  0.5 

— 

1  x  0 

-=-  /  q'dx,  F 

VIC  J  1  9 

°p  0 

31.4 

(computed ) 

— 

w,  lbs/hr 

— 

— 

+ 

0.02 

k,  Btu/(hrft2oF/f t) 

— 

— 

+ 

0.02 
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Run  No.  b,  Position  14" 


Variable 

Magnitude 

Unvertainty 

Interval 

Relative 

Uncertainty 

Gross  Power, 

Btu/ln 

27.51 

.  _  _ 

0.01 

Heat  Leak, 

Btu/in 

15.52 

— 

0.04 

(T„  -  T0),  °P 

1108 

5.0 

— 

/  q'dx,°F 
p  0 

759  ' 

(computed ) 

— 

w,  lts/hr 

— 

— 

0.02 

k,  Btu/(hrf t2oF/ft) 

— 

_ _ 

0.02 

Employing  the  root-mean-snuare  approximation  for  single 
sample  experiment  uncertainty  estimation,  the  following  prob¬ 
able  relative  uncertainty  was  computed  for  the  Nusselt  number 
determination . 

Run  No .  1  ±0.11 

Run  No.  6  ±0.11 

Thus  the  uncertainty  In  the  Nusselt  number  determination  is 
estimated  to  be  approximately  ±  11  per  cent  for  all  of  the 
test  runs. 

This  rather  large  experimental  uncertainty  is  in  major 
pa^t  attributable  to  the  large  relative  importance  of  the 
heat  leak,  and  the  ±  4  per  cent  uncertainty  placed  on  the 
heat  leak  calibration.  It  is  worth  noting  that  of  the  six 
test  runs,  in  only  two  do  the  measured  Nusselt  numbers  dif¬ 
fer  by  as  much  as  11  per  cent  from  the  theoretical  constant 
property  solution.  It  is  thus  felt  that  ±  11  per  cent  may 
be  quite  conservative.  Furthermore,  it  is  probable  that 
the  heat  leak  error  is  consistently  in  one  direction  rather 
than  random,  so  that  when  comparing  one  test  run  with  another 
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to  detect  effects  of  temperature  dependent  properties  the 
uncertainty  in  the  comparison  would  be  substantially  less 
than  +  11  per  cent. 


i 


Laminar  and  Turbulent  Cooling  Experiments  --  The  major 
sources  of  uncertainty  in  the  cooling  experiments  are  in  the 
measurement  of  the  inlet  and  exit  air  temperatures,  the  air 
flow  rates,  and  the  uncertainty  in  the  thermal  ’conductivity 
of  the  air.  The  difficulty  with  the  inlet  temperature  mea¬ 
surement  was  particularly  significant.  In  order  to  satisfy 
the  assumption  of  a  fully  established  velocity  profile  it 
was  necessary  to  have  a  velocity  developing  section  after 
the  inlet  air  temperature  measurement.  The  developing  sec¬ 
tion  was  wrapped  with  a  nichrome  ribbon  guard  heater,  and  a 
number  of  thermocouples  were  attached  along  the  tube  so  that 
the  guard  heater  could  be  adjusted  to  maintain  the  tube  wall 
temperature  the  same  as  the  inlet  air  temperature.  It  turned 
out  to  be  difficult  co  maintain  the  entire  tube  at  the  same 
temperature,  and  the  best  that  could  be  done  was  to  make  the 
temperature  versus  position  pattern  along  the  tube  approxi¬ 
mately  similar  for  each  run.  In. view  of  this  difficulty  it 
is  estimated  that  the  uncertainty  interval  in  temperature 
difference  between  the  inlet  air  and  the  cooling  water  is  no 
better  than  ±  5  per  cent.  With  this  temperature  difference 
reaching  as  high  as  1400°F,  the  uncertainty  interval  is  thus 
estimated  as  high  as  ±  70°F,  which  is  approximately  the  kind 
of  variation  that  was  measured  along  the  tube  for  the  high¬ 
est  temperature  runs.  The  uncertainty  in  the  exit  tempera¬ 
ture  is  considerably  less,  and  a  i  2  per  cent  estimate  in 
the  exit  temperature  difference  is  probably  quite  conserva¬ 
tive.  As  in  the  laminar  heating  experiments,  the  uncertainty 
in  flow  rate  and  thermal  conductivity  is  estimated  at  ±  2 
per  cent. 

With  these  estimated  uncertainties,  the  uncertainty  in 
the  resulting  Nusselt  numbers  are  calculated  as  follows  using 
the  root-mean-square  approximation. 
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±  0.048 
±  0.045 


l 


Laminar  Cooling  Experiments 
Run  No.  3 

Run  No.  13 

Turbulent  Cooling  Experiments 

Run  No.  2  ±  0.035 

Run  No.  23  +  0.064 

It  thus  is  concluded  that  for  the  laminar  cooling  ex¬ 
periments  the  uncertain cy  in  the  measured  Nusselt  numbers 
is  approximately  ±  5  per  cent,  and  for  the  turbulent  cooling 
experiments  the  uncertainty  is  approximately  +  6^  per  cent. 
These  results  are  consistent  with  the  scatter  of  the  experi¬ 
mental  points. 
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Ccff/r  Jvn#*rr  Con. 
vt  .  .err<.n#  Or, la 
(l)*ttr:  #r.  H,L.F«,y«r, 

Vie*  Ire*. 

Lorr-vij  Awrcra-tleel  l*i  . 

-w'l  %  >'*!*•  Si. 

»  ,  H  •  tL.1  • 

(i}At*  Elru.  T.l  •’*rs ,  l.  J 1  r-rle 
( 2  }  let  roll  £ll  icr.  Co. 

Sue  leer  Pow*r  Dl/Ia. 

Detroit,  .#:crixon 
Dr  »jr«r.H.4',eor!,  Ir^; . 

Boa  /.’*■,  Termini  Ante* 

Jo*  Argclee  Call'. 

(  1  )  At  t.r  j  Mr.  Gordon  M,J»c-sor, 
/Ice  PrealJ*r,t 

El tal-KcCul lough,  lr<  , 

"'*u  3.sr  M-i  t  *0  Ave. 
oar  Hruruj,  Calif, 
f  1  1  *  *  t  r  r  Mr.  J.S.VCul  lough, 
*yr.  Heaearrn 

Electric  Boat  Cornua ry 
<10 ton.  Connecticut 
(l)Atlni  Mr  .E.  S,  Derinlsr.n 

Elliott  Coepary 
Jtar.neltf,  renr.ay lvanl :i 
( l ) At tr,  t  Cr.  J.Ro.Jger  Shields, 
Dev. Engrg. Dept. 
PMrr.r.lld  Engine  and 
Airplane  Corp. 

Stratcs  DlvlaK-p 

Bay  Rhi»re,  Long  l*iard,N.Y. 

( 1 )  At  tn  t  hr,  V .  L.  a*hl  t  ney ,  ,’r . 

Rea.  Engineer 
p’errctfierw  Company 
1861  E,  65th  Stroirt 
Claveland  J,  Ohio 

(2) Attrx  Mi*.  Sven  Holm 

The  Pi«cr  Ccrp  .  Ltd, 

Research  and  Dev.  Dlv. 

P.O.Box  ^10 
Whittier,  Calif. 

(l)Attn:  Mr,  0 . 1. ,  Blohti, 

Mgr.  Re  see  re  h 
Ford  Hoto"  Coepany 
Lear borne,  Mlrnlger 
(l)Attrs  Hr.f?*ul  Kioto 7-h, 
h#e .  Cl<. 

( 1 )  Dj  ,C. N. (’rey , 

5c  lent  If  id  lat 

(1)  Mra.R/ieh*J  racbonall, 

Eng. Staff  LI  hi  cry 
;1)  Mr a.L.B..  hill Ipa, 

Sup&rlnt endept , 

Te,.h.  Ini *o. Sec . , 

Se  1  wr.fr  1  Me  lu*t 
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' '  '  i  •  J,  *»  ~  t  -  4p**- 
t  •  >l  *  i  *.*.  l.  Ux«# 

•  *#r*t*  for*, 

tr*»rjr\i»  fit  til** 

*  »*  i  rim  u*  :•  t, 

~*r  f:*»*ct*ia,  r*llf0f«|# 

>**r*l  V*.nr»  Corf  • 
c*  ttiut«r 
r.i*4f<lli  t ,  irattA-o 
.  )A*t«“t  **,  R.  ♦  . 

Clr,  et  lr«-g, 

■Jrrrral  *o!»r»  Ccrr. 

Rra#arch  Lid,  Mitator 
:«*irol t  J,  *ts»Usr. 

)A**ft  *.  .  #.  A,  fVrwrr.. 

4a  1  VtlrAi  f*pt. 
ti)  ir.  V.  C.  >*th 

<3er*-aJ  Motors  Ccrr. 
r^rrlron  bttiUr  i!i. 
Loc*p«rt,  S-r«  Yen 
l  •  )  A '  ’ ’r  ;  Mr .  J  uhr  • .  Jcslf  •;/ 

TV*  Jr  l  ACO-***u«»*  1  Co. 

.’*“  Mail  tor  A.'eri*. 

■«•■•  it>r*  */,  *«•  or* 

(*.  JAM,  ;  *,-.  S.  •.  » 

juir  nu  <*?-?. 

L - 1 '  ,  "in-/ 

ri!t,tur«rr  JO,  njylva-t* 

( 1 ) At ' '  t  Mr.  Charles  y.  *ott-a«v 
Tr<?  -*^jl -X-rrum^c  Co.,  lrc, 
Rrewster,  Ne«  Yo: 

(l)Att^t  Mr.  W.  .1,  *y  *>ov A r , 

/Ice  President 
Rubles  Aircraft  Company 
Florence  and  Traie  St*. 

Culver  City,  California 
(l)Attnj  Dor.a  H.  Conrsor,  Dir. 

Technlctl  Lltrary 
lr<*,Tcol  1 -Hanl  Cnr^mny 
11  Broadway 
K«rf  Yori<,  ?ieK  York 

«;!*'»»•?  »r.  *.  <*,  •viv^o*’ 

i rire r sol  1-Rand  Corf^ry 
n*  ;  1 1  lpsturg.  New  Jaraey 
(t'A'tr?  Dr.n.C.Dear,,  Jr, 

AJvancerl  ErgrR. 

(1)  .Kr.  J.  P.  Johnston, 

Conder.sor  and  Ejector 
Ergrtt.  Dept. 

International  Harveater  Co. 
Engineering  Reaearch  Dept, 

5?2$  D,  keatern  Blvd. 

Chicago  9,  Illinota 
(i)Attm  Dr.  Simon  K,  Chen 
I-T-E  Circuit  Breaker  . 
tOl  E.  Erie  Avenue 
Philadelphia  )U,  Pa 
(l)Attr.:  Mr.  Ashton  T.  Scott 

JjcU  and  Hetntz,  Inc. 
Engineering  Dept,,  Plant  4 
Cleveland  T,  Onlo 
(MAtint  Mra,  Harriet  D.  TriVillo, 

H.  W.  Kellogg  Company 
711  Thlr.l  Avenue 
New  York  17,  New  Vo*”.* 

(l)Attns  Mr.  P.  W,  Pvte'flon 
(l)  Mr,  0,  P.  laihvnbrenner 

(  1  Mr.  Ronald  E.  Smith, 

Vice  President 

Tl.e  Krai  as  L  Co.,  Inc. 

//J  WlllUmo  Avenun 
llackensack.  New  Jersey 
(l)Attni  Mr,  P.  Kralaal,  Jr. 

Lockheed  Aircraft  Co. 

Technical  Library 
Ikirtaifk,  “California 
(ijAttr;  p.  L.  .Meaelnger, 

Ei'grg.  No.  ?2?b 
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"  *•  *-+?••  f  *||  .M  ••  if*. 

•  M**»t 

CAtjAll, 

A**- 1  lr,  C.  >. 

•hill**  *a*  •? eit-wrt*#  f*. 
;V£  >!*-,• 

•'.Mlf*'.* 

UlAtT**!  Mr.C.  ?.  P*r*IfA#Pr*i  »  . 

V;r*e  *T«4-!t*n  i#J*  *  ••«  . 

Urj<r* 

M  Ah4*:*A  *1,,  C a  1 1 f , 

(lJAttrt  Mr.  •.  A.  JuUln 

fwrf*  i  Corfi.1  ■!  1  or 

Cv|»^  A,«r4* 

•5l»*«i*t  7,  4*.»roeA|n 
(IjAttrt  ^*.  •.  a.  Z-cf.l 
r l  #oi  farter 
•!J::as  street 
iW»  isf ■  Tor* 

i  Attn  M.*,C,T.*otlrtcr, 

Shell  Dwielcp-er*  CcrAjeny 
E»er/vlllA,  CAllfomla 
(l)Attrt  Dr.  C.  *.  Oerlet; 

A.  C.  2h!*>  Ccrp. 

ViJ'  Street 

)".l»ij/»e  J,  eUcor.jjr 
l  1  Uttnj  .  ».  h.  Mwr»r# 

Sole.-  Aircraft  Co*rar.y 
Sen  ytegj  U,  rell* 
f?)Attr i  Mr. T. A. fit! ,  CM»f 
Engineer 

(i)  Mr.  E. A. Crury.fi  oje 

F^glr«er,Dept .  eij 

Ctoleer  CMvelopner.t  Co 
>03  Wood tide  Avenue 
Issexvllle,  Michigan 
(1) Attn t  Hr,  t. A. Stalker 
Stanford  Reaearcn  lnat, 

Menlo  Terk,  California 
(IjAttn;  Dr.Nevlr  K.Kleiter 
1)  Mr. 5, H, Clark 

Stewart-M»rr<er  Corp. 

1  *S  1  Drover  Street 
lrkiianapolla  7,  Indiana 
(i)Attn:yr.R.D.Rardall , 

Mgr.  Research 

Sul zer  Bros,  Ltd 
*0  Church  Street 
Hew  York  ?,  New  York 
(l)Attn:  Mr.  Richard  Herold 
(i)Sverdrup  and  Percel,  Inc. 
Syndicate  Truat  Bldg. 

St.  Lcula  1,  Mlaaouri 
Texas  Eastern  TVanamlaalcn 
Corp. 

P.O.Box  1G12 
Shreveport,  Louisiana 
(l)Attn:  Hr.  C,  W.  Marvin 
Thompson  Products,  Inc. 
Acceaao“lea  Division 
2355  Euclid  Avenue 
Cleveland  17,  Ohio 

( 1 )  Attn t  Mr. 3. B. Spangler, Jr. 
Tho-npson  Products,  Inc, 
Research  and  Development 
Jet  Dlvlalon 
Willoughby,  Ohio 

{IjAttn*  Mr. W. A.Corpton, Krg . 

The  Trane  Company 
2nd  and  Cameron  Avenues 
•  La  Croase,  Wisconsin 

(2) Attn :  Mr.  H.  C.  Rooka 

(1)  Mr.  I.  T.  Wetzel 

(?)Unlon  Carbide  Nuclear  Co. 
hak  pidge  Gaseous  Diffusion 
Plant 

Plant  Records  Department 
P.  0.  Box  P 
Oik  Ridge,  Tennessee 
United  Aircraft  Corp. 

^00  Hein  Street 
Ea3t  Hartford  8,  Conn, 

(l)Attns  Hr.  Holvrt  C.  Sale, 
Chief  Librarian 
(1)  Hr.  Ned  C.  Ric*r  Jr., 

Research  Dept. 
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MiAirtt 

Cl)  U»ll  iroc.eor 

I  >*ri  *j  1 1  a t or  Co«^*-y 
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WttCOTAln 

(i)Af-t  r-,  H-.f.Prir^r, 

Asa.  engineer 

Yu  .»  VA«arcn  •»>i  >,ei. 
i  ertsr 

VtrlA**  «•, 

Ui  Carlo*.  Caivr. 

(1  A*trf  A*,  tjf  rT* ! r  £.  lia.-xl# 

(  i  )Hr .  IW;*or  k.  4l.*r, 

.*>>*  All  jew  JI^Ai 
Her.l*  f*P4,  CaMr 

( 1 ) Mr ,  i.  R.  Bgtrook 
Vic*  F resident 
Nucor  Reirtrch,  Iiw;, 
ie^il  eolcott  Avenge 
Mrrnd lie  ?C,  Mich. 

U)Mr.  P.  A.  Brocks 
r>pr,  if  Agr.  Er.grg. 
tfr.lv.  of  Ctllf.  Agr. 

Exi-ers  vr  t  314110*“ 

Davia,  Caltforr.!* 

( 1 ) Mr .  w.  I.  Comtole 
?9 -Brentwood  Drive 
Bloomfield,  Correct icut 

MlHr.  Paul  Dawson 
?^«ro  cu>» .  Sallrlath 
Clnclnra'  1  UJ,  Ohio 

(l)Mr^  A.  J .  Ede 

Heat  DlvUlon 
Me  »h.  Erg.  Reaearch  Lot 
D.l  . I.R. ,  Kilbride 

neai  Glsstow,  SCOTLAND 
SEHI  IDl 

Office  of  the  Asat.  Naval 
Attache  for  Research 
Naval  Attache,  Arerlcan 
Embassy 

Navy  No.  100,  Flaat  P.  0. 
New  York,  New  York 

( 1 ) Hr .  C.  M.  Flxman 

C.  tier  Marine  Engineer 
Code  25^ 

Ran  Pranclaco  NavbI 
Shipyard 

San  Pranclso  25»  Calif, 
(l)Mr.  Ralph  M.  Hclntz 
P.  0.  Box  5*6 
Los  Catoa,  Ca.'if. 

( 1 J  Mr .  Altart  L.  Holiday 
550  Lakewood  Clrcla 
Walnut  Creek,  Itall'. 

(l)Mr.  Robert  V.  Hlelnachmldt 
20  East  Street 
Stone  ha  as  80,  Ka.is, 

(l)Rear  Adm,  Paul  P.  Lee, 

U3N  Ret. 

1^975  Manuel la  Avenua 
Loa  Altoa,  Calif. 

(1 )Mr.  E.  J.  Le  Fevre 

Department  of  Mechanical 
Engineering 
Queen  Mary  Collage 
London  E.  1.  ENGLAND 
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